Stability of the superfluid state in a disordered ID ultracold fermionic gas 



Masaki TezukeQ 

Department of Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japai^\ 



o 

(N 
-t— > 
O 

O 

(N 



<Z3 

^— > 

•4-* 



c 

o 
o 



> 
m 

(N 
(N 

(N 

o 



X 



Antonio M. Garcia-Garci'a 
CFIF, Instituto Superior Tecnico, Universidade Tecnica de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 

We study a ID Fermi gas with attractive short range-interactions in a disordered potential by the 
density matrix renormalization group (DMRG) technique. Our results can be tested experimentally 
by using cold atom techniques. We identify a region of parameters for which disorder enhances 
the superfluid state. As disorder is further increased, global superfluidity eventually breaks down. 
However this transition seems to occur before the transition to the insulator state takes place. This 
suggests the existence of an intermediate metallic 'pseudogap' phase characterized by strong pairing 
but no quasi long-range-order. 

PACS numbers: 67.85.Lm, 67.25. dj, 37.10.Jk, 72.15.Rn 



I. INTRODUCTION 

It is now possible to realize experimentally disorder 
and interactions with unprecedented precision by using 
cold atom techniques 0, 0. This is an ideal setting 
to test theoretical predictions on novel phases of quan- 
tum matter and quantum phase transitions Q. Moti- 
vated by these possibilities we study a disordered ID 
Fermi gas with short-range attractive interactions by the 
DMRG technique. The effect of disorder is mimicked 
by a quasiperiodic (multichromatic) potential. Both the 
potential [2] and the interaction can be implemented ex- 
perimentally. Our main results can be summarized as 
follows: a) attractive interactions enhance localization ef- 
fects. The critical disorder at which the metal-insulator 
transition occurs decreases as the interaction becomes 
stronger; b) in contrast to higher dimensions, fluctuations 
in the metallic phase, but close to the insulator transi- 
tion, break down quasi long-range order. The resulting 
anomalous metallic region has 'pseudo-gap' features; c) 
in the superfluid phase, and for moderate interactions, 
disorder enhances quasi long-range order. 
We start with a brief overview of previous research on 
this problem. In the non-interacting limit the nature of 
the eigenstates of a ID tight-binding model, with hop- 
ping t = 1, in the quasiperiodic potential [2| 



V(n) = Acos(27ra;n + 6) 



(1) 



with oj irrational and £ [0, 2tt) depends on the value of 
the disorder strength A > 0. All the eigenstates are ex- 
ponentially localized [1, Q for A > 2 with a localization 
length oc 1/|A — 2|. For A < 2 the quantum dynamics 
is similar to that of a free particle in a periodic poten- 
tial. For A = 2 the system undergoes a metal-insulator 
transition [5]. We note that the potential is strongly 
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correlated (V(n)V(0)) oc cos(27rwn) 7]. In ID, a non- 
decaying (V(n)V(0)) is a necessary condition |6| for the 
existence of a band of metallic states. In the limit A — >• 
an exact solution for a continuous ID model with short 
range attractive interactions - the Gaudin - Yang model 
I8J - is available [8 10]. An exact solution is also known 
for the discrete version of this model, the ID Hubbard 
model PQ. For \U\ <C 1 pairing is BCS-like. For 
— > oo the system behaves as a hard-core Bose gas 

It was found in [14| that the addition of a weak Gaus- 
sian disorder induces a metal-insulator transition for suf- 
ficiently strong interactions. The effect of a quasiperiodic 
potential has also been addressed in the literature flEh 
flSj ] . The numerical results of [15J indicate that the critical 
disorder at which the metal-insulator transition occurs 
depends on the strength of the interaction. By contrast 
the DMRG analysis of [n| concluded that, for spinless 
fermions, the critical disorder is the same as in the non- 
interacting case. In [ijj], also employing a DMRG tech- 
nique, it was found that the presence of a weak disordered 
potential enhances superfluidity 

Bosonization techniques combined with a renormaliza- 
tion group analysis were employed in [18J to investigate 
the effect of interactions in another ID quasiperiodic sys- 
tem, the Fibonacci chain [2(J. The perturbative treat- 
ment of [HI showed that the critical disorder depends 
on both the strength of the interactions and the position 
of the Fermi level. We note that, as in Eq.([T]), corre- 
lations of the potential studied in [lj| are very strong 
lim n - i . O o(V(n)V(0)) ^ 0. However different reasons pre- 
vent from a direct comparison between these models: a) 
in [IH only spinless fermions are considered, b) for A <C 1, 
the limit in which the formalism of [l8| is applicable, our 
system is in the metallic region with properties almost 
identical as those of a periodic potential, c) for no in- 
teractions the spectrum of the Fibonacci chain is singu- 
lar continuous for all A. This leads to eigenstates that 
are power-law localized and quantum superdiffusion [2(| . 
Such features are only found in Eq.(fT]) for A = 2. For 
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A < 2 the spectrum of Eq.([T]) is absolutely continuous as 
in a perfect metal. 

For results on the dynamics of a Bose gas in a 
quasiperiodic potential we refer to [2l|. Mean field ap- 
proaches in ID are problematic since fluctuations, spe- 
cially in the presence of a disordered potential, are not 
negligible. We thus anticipate qualitative differences with 
respect to the 2D and 3D cases where, for disorder weak 
enough, quasi long-range order persists [22| even in the 
insulator region provided that the localization length is 
larger than the coherence length. Finally we mention 
that the effect of disorder in Fermi gases of higher di- 
mensions has been investigated in [231 ] using mean field 
techniques and neglecting Anderson localization effects 
[2H . For numerical studies on the attractive Hubbard 
model in a disordered potential we refer to [251 ] . 



II. THE MODEL 



We study the discrete L-site Hubbard model, 



L-l 



L-l 



H = - ^ (c|_ li0 .Ci, ff + h.c.) + U ^ ^,t"»4 

i— 1.(7 2—0 

L-l 

+ J2 y (*)"*> 



(2) 



i=0 



where Ci y(T annihilates an atom at site i in spin state <r(=t 

,4), n i)(T = cj,oA^ = r \t + ™a> v(i) is s iven by 

Eq.(TT]) with uj = F n -i/F n the ratio of two consecutive 
Fibonacci numbers, L = F n + 1 and 9 = so that V(0) = 
V(L — 1) = A. We note that we have set the hopping 
integral t = 1. 

The behavior of Hamiltonian Eq.@ in certain limits 
is already known: a) for \U\ 3> 1 the system maps onto a 
weakly interacting bosonic gas with a kinetic term which 
is 1/\U\ smaller than in the original fermionic model. 
Therefore the critical disorder at which the transition 
to localization occurs is A c « 2/|c7 26]; b) the coher- 
ence length for weak disorder (A <C A c ) is £ co oc 1/U 2 for 
\U\ > 1 and Ceo oc e 1 /!^! for \U\ < 1; c) for U < 1 not 
very large, the spin gap (see Eq.©) As oc l/£i oc , HI] 
with Ci oc the localization length. 

The above information is enough to put forward a ten- 
tative description of the system phase diagram (in the 
U < 0, A plane): a) for fixed \U\ > 1 and \U\ < 1, the 
loss of quasi long-range order and the transition to the 
insulator phase will occur at similar A's: A c ~ 2/\U\ and 
A c ~ 2 respectively; b) for intermediate U it might be 
possible that the two transitions take place at slightly 
different A's as the breaking of superfluidity might be in- 
duced by phase and amplitude fluctuations in the metal- 
lic region. 

In order to test the validity of these qualitative argu- 
ments we study the Hamiltonian Eq.([2]) with the poten- 
tial Eq.© by the DMRG technique. The filling fac- 
tor v is kept constant v — N/L = 1/9 for (N,L) — 
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FIG. 1. (Color online) Upper: Pairing structure factor P s 
Eq.Q, as a function of the system size L for different A's 
and U = —6. Superfluidity, characterized by an increasing 
P S (L), is observed up to A c ~ 0.29. Lower: P s as a function 
of disorder also for U = — 6 and different sizes. A P s almost 
independent of L is a signature of broken quasi long-range 
order. See text for more details. 
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FIG. 2. (Color online) Same as Fig[T]but for U = -1. Upper: 
P s only increases with L for A < 0.95. Global superfluidity is 
thus broken at A c w 0.95. Lower: P s is an increasing function 
of A until A ss 0.8. Therefore the quasiperiodic potential 
enhances superfluidity for moderate disorder. 



(10,90), (26,234), (42,378) - quantitatively our results 
might depend on v <C 1 [18| -. We obtain the ground 
state for N spin-up and N spin-down atoms. Up to 
to = 400 basis states for each block are kept in the finite- 
size system DMRG iterations. 
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FIG. 3. (Color online) Upper: Ie, Eq.Q as a function of L, 
for different A's and U = —6. A metal-insulator transition 
is observed at Aj" s « 0.31. However for [/ = 0, A^ ns = 2. 
Therefore attractive interactions enhance localization. Lower: 
Ie as a function of A. An increase of Ie with the system size 
is a signature of a metal. 
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FIG. 4. (Color online) Same as FigEJbut for U = -1. Upper: 
The metallic state is characterized by a Ie that increases with 
L. The insulator transition occurs at AJ. 118 ~ 1.0. In contrast 
to the U = —6 case, it is observed a further increase of Ie 
very close to the transition A m 0.99. This is a consequence of 
the enhanced eigenfunction correlations in this region 0, [2S| . 
For U = — 6 the coherence length is much smaller and con- 
sequently eigenfunctions correlations are suppressed. Lower: 
Also, in contrast with the U = —6 case, the metallic state is 
also enhanced for intermediate disorder A's below the transi- 
tion. This is again a quantum coherence effect caused by the 
bands structure of the quasiperiodic potential (see text). 




FIG. 5. (Color online) Spin gap As, Eq.© as a function of 
A for different U's and fixed L. For small A the gap is an 
increasing function of disorder as a consequence of the band 
structure of the quasiperiodic potential. Close to the insulator 
transition A < A c there is an additional gap enhancement 
caused by eigenfunction correlations 



III. RESULTS 

Our first task is to determine for what range of pa- 
rameters global superfluidity breaks down. In weakly 
disordered BCS superconductors a study of the ground 
state and the low energy excitations is enough to answer 
this question as the vanishing of the gap is equivalent to 
the breaking of global coherence. In strongly disordered 
and strongly coupled superconductor the situation is dif- 
ferent as the gap might be finite even after fluctuations 
have destroyed global superfluidity [22| ■ It is thus nec- 
essary to compute observables that directly measure the 
phase stiffness of the system. 



A. Phase rigidity: Pairing structure factor 

A popular choice [25| is the averaged equal-time pairing 
structure factor, 



(3) 



where (...) stands for the spatial average on the site in- 
dex i, T(i,r) = (A(i + r)A' l '(i)) and A(j) = CjfCj^. Quasi 
long-range order (there is no true order in ID) occurs for 
r(r) ~ l/r K for r 3> 1. In the case with no disorder it 
was demonstrated in [l2| that superconductivity corre- 
lations are always leading with respect to other types of 
quantum order and that K < 1. The limit K — 1 is only 
achieved in the limit \U\ — > oo. In the disordered case it 
is also plausible to expect that K < 1 but K = 1 occurs 
for a finite U which depends on A. Therefore we define 
quasi global superfluidity by P s cx L 1 ~ K with K < 1. In 
Figs. [T]and[2]we observe: 

a) the critical A = A c < 2 at which global superfluid- 
ity breaks down decreases as \U\ increases. Therefore a 
tighter binding is correlated with a greater instability to 
disorder effects [26| . 

b) For not too strong U, P s is an increasing function of 
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A up to some A close but smaller than A c . 
c) For \U\ S> 1 this feature is not observed. 
We believe that b) is a coherent effect related to the pe- 
culiar band structure induced by the quasiperiodic po- 
tential. This is also consistent with c). As \U\ increases 
the coherence length decreases, the details of the spec- 
tral density are smoothed out, and no enhancement of 
superfluidity is observed. 



range of A's for which we observe this phase is relatively 
narrow and it seems to decrease for larger U. Therefore 
we cannot discard the possibility that this metallic phase 
is a finite size effect, namely, the system is already an 
insulator but the localization length is larger than the 
system size. 

C. Low energy excitations: Spin gap 



B. Localization: Density fluctuations 

We now turn to localization properties. More specifi- 
cally we determine numerically the location of the critical 
disorder A" ls at which the metal-insulator transition oc- 
curs. Different quantities, such as density fluctuations 
[26j or the conductance [l6j], provide a similar estima- 
tion of localization effects. However the numerical value 
of AJ. ns might depend weakly on the observable employed 
[27j . We present results for the density fluctuations, 



Finally we study the low energy excitations of Eq.© 
by computing the minimum energy to break a pair, the 
so-called spin gap, 



I E = \Y, Sn (hN,N) 2 



(4) 



where Sn(i, N, N) = n(i, N + 1, iV + 1) - n(i, N, N) is 
the ground-state atomic density at site i for N spin-up 
and N spin-down atoms, E stands for the ground state 
energy in this case. For U — 0, it corresponds with the 
usual definition of the inverse participation ratio in non- 
interacting systems (28j. 

In the insulator region it is proportional to the localiza- 
tion length Ie oc £i oc . It decreases slowly as disorder 
increases until it saturates Ie — > 1/4 for A — > 00. In 
the metallic region (A <C AJ. ns ), Ie oc L with only a weak 
dependence on A. Close to the critical region, Ie oc L a 
with a < 1 a constant that depends on the eigenstates 
multifractal dimensions [28| . 

In Figs. [3] and 0] it is shown that A™ s decreases with 
\U\. This enhancement of localization effects caused by 
attractive interaction is consistent with previous results 
in the literature (2(| . It is also observed that for U = — 1 
the dependence on A is not monotonous. Initially it de- 
creases with A but close to the transition (A < AJ. ns ) has 
a sharp peak before a steep drop right at A™ s . This is 
not expected as it is believed that quasi long -range or- 
der is always weakened by disorder effects [22| ■ Within a 
mean field approach this might be attributed to the en- 
hancement of eigenstate fluctuations around the critical 
region Q . The absence of enhancement for larger \U\ is a 
consequence of the shorter coherence length in this case. 
Single particle fluctuations are suppressed if the coher- 
ence length becomes smaller than the system size. 
We note that, according to Fig. [H the transition to lo- 
calization occurs at AJ. ns ~ 1.0. On the other hand, ac- 
cording to Figj2] , global superfluidity breaks down at 
A c ~ 0.95. This suggests the existence of a metallic 
pseudo-gap phase for 0.95 < A < 1.0 characterized by 
strong pairing but no global superfluidity. We note the 



A s ee E (N + 1, N - 1) - E (N, N), 



(5) 



where Eo(Nf, N±) is the ground state energy for Nf spin- 
up and iVj, spin-down atoms. In Fig. [5] we present re- 
sults for As for a fixed I as a function of A and different 
IPs. It is observed that As is an increasing function of 
A. By contrast in 2D weakly disordered systems the gap 
decreases with A 22[ since the spectral density around 
the Fermi energy decreases with disorder. In quasiperi- 
odic systems the situation is different. As A increases, 
the spectral density around the Fermi energy develops 
gaps at different scales and the spectral density in the 
remaining bands becomes higher. For not too large A's 
it is likely that, on average, there are no gaps around the 
Fermi energy. Therefore both the spectral density and 
the spin gap will increase as A increases. Close to the 
metal-insulator transition, strong density-density fluctu- 
ations in the one-body problem 0, further enhance 
the gap. This enhancement is a coherent effect and there- 
fore it is expected to diminish as the coherence length 
becomes of the order of the system size which occurs in 
the region of strong coupling. For larger A, already in the 
insulator region, the spin gap As still increases with A. 
This is not related to superconductivity but rather to the 
fact that now the gap is related to the mean level spacing 
which in the insulating region increases with disorder Q . 



IV. CONCLUSIONS 

We have studied the stability of the superfluid state 
in a ID interacting and disordered Fermi gas. We 
have shown attractive interactions enhance localization 
effects. For intermediate couplings \U\ « 1 we have 
identified a region close to the insulator transition in 
which superfluidity is substantially enhanced. Moreover 
our numerical results suggest that the breaking of global 
superfluidity might occur at a slightly weaker disorder 
than the insulator transition. If this is confirmed, a 
"pseudo-gap" metallic region characterized by pairing 
but no global superfluidity occurs between the two tran- 
sitions. These results provide a theoretical framework 
for experimental studies of quantum phase transitions in 
ID cold Fermi gases. 
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